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Abstract 

In order to identify the most favorable situation for superconductivity in 
the repulsive single-band Hubbard model, we have studied instabilities for 
d-w&ve pairing mediated by antiferromagnetic spin fluctuations and p-pairing 
mediated by ferromagnetic fluctuations with the fluctuation exchange approx- 
imation in both two dimensions and three dimensions. By systematically 
varying the band filling and band structure we have shown that (i) ti-pairing 
is stronger in two dimensions than in three dimensions, and (ii) p-pairing is 
much weaker than the (i-pairing. 
PACS numbers: 71.10.Fd, 74.20.Mn 
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The discovery of the high-temperature superconductivity in copper oxides by Bednorz 
and Mullen! has kicked off intensive studies for electron mechanisms of superconductiv- 
ity. Specifically, it is becoming increasingly clear that superconductivity can arise from 
repulsive electron-electron interactions. A persuasive scenario is that the superconductivity 
comes from a pairing interaction mediated by antiferromagnetic(AF) spin fluctuations. A 
phenomenological calculationlHl along this line has succeeded in reproducing anisotropic d- 
wave superconductivity as well as anomalous normal-state properties. Analytic calculations 
on a microscopic level with the fluctuation exchange approximation (FLEX), developed by 
Bickers et a/.i, has also been applied to the Hubbard model on the two-dimensional(2D) 
square lattice!'! to show the occurrence of the superconductivity. Numerically, a quantum 
Monte Carlo study has indicated the pairing instability!!. 

These results indicate that the superconductivity near the AF instability in 2D has 
a 'low Tc ~ O(0. Oli) it: transfer integral), i.e., two orders of magnitude smaller than 
the original electronic energy, but still 'high Tq ~ 0(100 K) for t ~ 0(1 eV). Then the 
next fundamental questions, which we address in this paper, are: (i) Is 2D system more 
favorable for spin- fluctuation mediated superconductivity than in three dimensions(3D)? 
(ii) Can other pairing, such as a triplet p-pairing in the presence of ferromagnetic spin 
fluctuations, become competitive? We take the single-band, repulsive Hubbard model as a 
simplest possible model, and look into the pairing with the FLEX method both in 2D and 
3D. The FLEX method has an advantage that systems having large spin fluctuations can 
be handled. 

Let us touch a little more upon the background to the above two questions. The pos- 
sibility of triplet pairing mediated by ferromagnetic fluctuations has been investigated for 
super fluid a heavy fermion system UPt 3 0, and most recently, an oxide Sr 2 RuO 4 It 

was shown that ferromagnetic fluctuations favor triplet pairing first by Layzer and Fay0 be- 
fore the experimental observation of p-wave pairing in 3 He. For the electron gas model, Fay 
and LayzerEl or later ChubukovS has extended the Kohn-Luttinger theoremS to p-pairing 
for 2D and 3D electron eas in the dilute limit. TakadaS discussed the possibility of p-wave 



superconductivity in the dilute electron gas with the Kukkonen-Overhauser mo deM As for 
lattice systems, 2D Hubbard model with large enough next-nearest-neighbor hopping (t') 
has been shown to exhibit p-pairing for small band fillingsJiil HlubinaS reached a similar 
conclusion by evaluating the superconducting vertex in a perturbative way.0 However, the 
energy scale of the p-pairing in the Hubbard model, i.e., Tc, has not been evaluated so far. 

As for 3D systems, Scalapino et showed for the Hubbard model that paramagnon 
exchange near a spin-density wave instability gives rise to a strong singlet c?-wave pairing 
interaction, but Tc was not discussed there. Nakamura et all extended Moriya's spin 
fluctuation theory of superconductivity! to 3D systems, and concluded that Tc is similar 
between the 2D and 3D cases provided that common parameter values (scaled by the band 
width) are taken. However, the parameters there are phenomelogical ones, so we wish to 
see whether the result remains valid for microscopic models. 

Here we shall show that (i) <i-wave instability mediated by AF spin fluctuation in 2D 
square lattice is much stronger than those in 3D, while (ii) p-wave instability mediated by 
ferromagnetic spin fluctuations in 2D are much weaker than the d-instability. These results, 
which cannot be predicted a priori, suggest that for the Hubbard model the 'best' situation 
for the pairing instability is the 2D case with dominant AF fluctuations. 

We consider the single-band Hubbard model with the transfer energy tij = t(= 1 here- 
after) for nearest neighbors along with = t' for second-nearest neighbors, which is included 
to incorporate the band structure dependence. The FLEX starts from a set of skeleton di- 
agrams for the Luttinger-Ward functional to generate a (/c-dependent) self energy based 
on the idea of Baym and Kadanonfi Hence the FLEX approximation is a self-consistent 
perturbation approximation with respect to on-site interaction U. 

To obtain Tc, we solve, with the power method^, the eigenvalue (Eliashberg) equation, 

AS(2) W = ^E S(2) (^)|G(A;')|V 2 )(A;-A;'), (1) 

where 
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for spin triplet pairing, where Xo(q) = —T/NJ2kG(k)G(k + q) is the irreducible susceptibil- 
ity, G{k) the dressed Green's function, and E( 2 )(fc) the anomalous self energy. At T = Tq, 
the maximum eigenvalue Am sx reaches unity. We take N = 64 2 sites with n c = 2048 Mat- 
subara frequencies for 2D, or iV = 32 3 with n c = 1024 for 3D. 

Let us start with the 2D case having strong AF fluctuations. In Fig.[l], we plot Xrpa(<z) = 
Xo/(l~ Uxo) as a function of the momentum for the 2D Hubbard model with t' — 0, n — 0.85 
(nearly half-filled) with U = 4 and T = 0.03. A dominant AF spin fluctuation is seen from 
Xrpa peaked near (vr,7r). 

We can then solve the Eliashberg equation ([J) to plot in Fig.|3|(a) Am &x as a function 
of temperature T (normalized by t). The behavior of \G(k, iirksT)^ that appear in the 
Eliashberg equation is indicated in Fig.|I|. How AMax is close to unity measures the pairing, 
and AMax tends to unity at T ~ 0.02, in accord with previous resultsi'i!. We also plot the 
reciprocal of the peak value of XRPA(k, 0), where 1/x — > indicates the magnetic ordering. 
While we cannot compare Ap^ax and Xrpa on an equal footing, since pairing fluctuations are 
neglected in the Eliashberg equation while the susceptibility is treated beyond the mean 
field, we can discuss the behavior of AMax when the situation is varied. 

Keeping the above result in mind as a reference, we move on to the case with ferro- 
magnetic spin fluctuations, where triplet pairing is expected. This situation can be realized 
for relatively large 0.5) and electron density away from half-filling in the 2D Hubbard 
model. Physically, the van Hove singularity shifts toward the band bottom with t', and 
the large density of states at the Fermi level for the dilute case favors the ferromagnetism. 
It has in fact been shown from quantum Monte Carlo study that the ground state is fully 
spin-polarized at t' = 0.47, n ~ 0.4.00 

We have calculated AMax for the density varied over 0.2 < n < 0.6 and t' varied over 
0.3 < f < 0.6 for U = 4, 6 with T = 0.03, and have found that Am &x becomes largest for 



n = 0.3, t' = 0.5, so we concentrate on this parameter set hereafter. If we look at in Fig. 
|2] the momentum dependence of \G(k, mk,BT)\ 2 and Xrpa for this case with U — 4, Xrpa is 
indeed peaked at Y (k = (0, 0)). The question then is the behavior of Am sx as a function of 
T, Fig.[5](b), which shows that Am &x is much smaller than that in the AF case, Fig.^a). 

A low Tq for the ferromagnetic case contrasts with a naive expectation from the BCS 
picture, in which the Fermi level located around a peak in the density of states favors 
superconductivity. We may trace back two- fold reasons why this does not apply. First, if we 
look at the dominant (oc 1/[1 — Uxo(q)]) term of the pairing potential itself in eqs. @ 
and (H), the triplet pairing interaction is only one-third of that for singlet pairing. Second, 
the factor \G\ 2 for the ferromagnetic case (FigfJ) is smaller than that in the AF case (Fig.^J), 
which implies that the self-energy correction is larger in the former. Larger self-energy 
correction (smaller |G| 2 ) leads to smaller eigenvalues of the Eliashberg equation ([!]). Even 
when we take a larger repulsion U to increase the triplet pairing attraction (susceptibility), 
this makes the self-energy correction even stronger, resulting in only a small change in A. 

Let us now move on to the case of <i-wave pairing in the 3D Hubbard model. In this case, 
we find that the representation of groupil has the largest AMax, so we look at this 
pairing symmetry hereafter. We have calculated AMax for the density varied over 0.75 < n < 
0.9 and t' varied over -0.5 <f < +0.4 for U = 4,6, 8, 10, 12 with T = 0.03. Among these 
parameter sets, we have found that AMax becomes largest for n = 0.8, t' = —0.2 ~ —0.3 and 
U = 8 ~ 10, so hereafter we concentrate on this parameter set. 

In Fig. 0(c), we again plot AMax along with the reciprocal of the peak value of XRPA(k, 0) 
as a function of T for t' = —0.2, —0.3 ,{7 = 8 and n = 0.8. We can immediately see that 
the pairing tendency in 3D is much weaker than that in 2D. Technically, for the sample size 
/V = 32 3 and the number of Matsubara frequencies n c = 1024 there are some finite-size 
effects for T < 0.02. As the inset for a larger n c = 2048 exemplifies, however, AMax tends to 
increase with N and n c , and we believe that a finite Tc (< 0.01) may be obtained at least 
for t' = —0.3, U = 8, n = 0.8 in the limit of large iV and n c , but this is still significantly 
smaller than in 2D. 
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Having confirmed this, the question now is: why is the ^-superconductivity much stronger 
in 2D than in 3D? We can pinpoint the origin by looking at the various factors involved in 
the Eliashberg equation. Namely we question the height of and \G\ 2 along with the 
width of the region, both in the momentum sector and in the frequency sector, over which 
V^ 2 \k) contributes to the summation over k = (k, iu n ). 

We first plot \G\ 2 for k z = 0, tt/2, tt as a function of k x and k y in the 3D Hubbard model 
for t' = —0.2, n = 0.8 with U = 8 in Fig. H We can see that the maximum of |G| 2 in 3D, if 
multiplied by U 2 arising in the Eliashberg equation, is in fact larger than in 2D. Were this 
factor the origin, a larger Amrx would result in 3D. 

We can then question how the peak in Xrpa spreads in the frequency axis. Fig. [|(a) 
displays Im;\;RPA(&Max, ^) (kMax : the momentum for which x(k, 0) is maximum) as a function 
of u (obtained by an analytic continuation with Pade approximation^!) . The figure compares 
the 'best 3D' case (f = -0.2, n = 0.8, U = 8) with a typical 2D case with t' = 0, n = 0.85 
and U = 4 having a similar magnitude of \. We can see that Inr\;(co>), when this quantity is 
normalized by its maximum value while u by t, exhibit surprisingly similar behaviors for 2D 
and 3D. So we can exclude the frequency width from the reason for the 2D-3D difference. 
Note that if the frequency spread of the susceptibility scaled not with t but with the band 
width, as Nakamura et a@ have assumed, Aiviax would have become larger. So this is one 
reason why we stress that the present result that 2D is the best is by no means readily 
predictable. 

If we turn to the momentum sector, Fig. [|(b) for XRPA(k, 0) shows that the width, a, of 
the XRPA(k, 0) peak in each momentum direction is similar to those in 2D (Fig.|l|). Since the 
right-hand side of the Eliashberg equation ([I]) is normalized by iV oc L D with L being the 
linear dimension of the system, A oc (a/L) D is smaller in 3D than that in 2D when the main 
contribution of to A is confined around (tt, tt) or (tt, tt, tt). So we can conclude that this 
is the main reason why 2D differs from 3D. 

We have also obtained results (not shown here) in 3D for the body centered cubic lattice 
near half- filling (where strong AF fluctuations are expected), but the <i-pairing is again weak. 
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The p-pairing in the face centered cubic lattice with low band filling (where ferromagnetic 
fluctuations are expected) is found to be even weaker. These results will be published 
elsewhere. 

To summarize, <i-pairing in 2D is the best situation for the repulsion originated (i.e., spin 
fluctuation mediated) superconductivity in the Hubbard model. In this sense, the layer-type 
cuprates do seem to hit upon the right situation. However, our conclusion has been obtained 
for the simplest possible single-band Hubbard model, while the detailed behavior of Tq may 
depend on the model. Indeed, if we turn to other 3D superconductors, the heavy fermion 
system, in which the pairing is thought to be meditated by spin fluctuations, the Tc, when 
normalized by the band width W, is known to be of the order of 0.0011U. Since the present 
result indicates that T c , normalized by W, is ~ 0.0001IU at best in the 3D Hubbard model, 
we may envisage that the heavy fermion system is an instance in which larger frequency 
and/or momentum spreads in x(k, oj) are utilized than in the Hubbard model. 

After completion of this study, we came to know the work by Monthoux and Lonzarich.il 
Using a phenomenological approach, they conclude for 2D systems that the c?-wave pairing 
is much stronger than p-wave pairing, which is consistent with the present result. 

We would like to thank K. Ueda and H. Kontani for illuminating discussions. R.A. 
would like to thank S. Koikegami for discussions on the FLEX. R.A. is supported by a 
JSPS Research Fellowship for Young Scientists, while K.K. acknowledges a Grant-in-Aid for 
Scientific Research from the Ministry of Education of Japan. Numerical calculations were 
performed at the Supercomputer Center, ISSP, University of Tokyo. 
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FIGURES 




(0,0) (Jt,0) (7t,7i) (0,0) 

FIG. 1. The squared absolute value of Green's function for the smallest Matsubara frequency, 
iuj n = mk^T (left) and the RPA spin susceptibility (right) against the wave number for the 2D 
Hubbard model with t' = 0, n = 0.85 and [7 = 4. 
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FIG. 2. A similar plot as in Fig. |l| for the 2D Hubbard model for a finite t' 
smaller n = 0.3 with U = 4. 
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FIG. 3. The maximum eigenvalue of the Eliashberg equation (solid lines) and the reciprocal 
of the peak of xrpa (either Ferro- or Antiferro-magnetic, dashed lines) against temperature for the 
Hubbard model in (a)2D with if = 0, n = 0.85 and U = 4, (b)2D with t' = 0.5, n = 0.3 and {7 = 4, 
(c)3D with t' = —0.2, —0.3, n = 0.8 and U = 8. The inset in (c) is the results for a larger number 
of Matsubara frequencies (=2048) for t' = —0.3. 
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FIG. 4. (a)IniXRPA(kMax> to) (normalized by its maximum value) as a function of tu/t for 
3D Hubbard model with t' = -0.2, n = 0.8, U = 8 and T = 0.03, 0.04 (dashed line) and for 2D 
Hubbard model with t' = 0, n = 0.85, U = 4 and T = 0.03, 0.04 (solid line). For T = 0.03 2D and 
3D results almost overlap with each other. (b)RPA spin susceptibility XRPA(k, 0) as a function of 
the wave number for 3D Hubbard model with t' = —0.2, n = 0.8, T = 0.03 and U = 8. 
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FIG. 5. A plot for Green's function against k x and k y with k z 
model with t' = -0.2, n = 0.8 U = 8, T = 0.03. 
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